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Abstract 

We present an algorithm for maintaining maximal matching in a graph under addition and deletion 
of edges. Our data structure is randomized that takes (9(log?i) expected amortized time for each edge 
update where n is the number of vertices in the graph. While there is a trivial 0(n) algorithm for edge 
update, the previous best known result for this problem was due to Ivkovic and Llyod|4]. For a graph 
with n vertices and m edges, they give an 0((n + m) ' 7072 ) update time algorithm which is sublinear 
only for a sparse graph. 

For the related problem of maximum matching, Onak and Rubinfeld J5) designed a randomized 
data structure that achieves (3(log 2 n) expected amortized time for each update for maintaining a c- 
approximate maximum matching for some large constant c. In contrast, we can maintain a factor two 
approximate maximum matching in O(logn) expected amortized time per update as a direct corollary 
of the maximal matching scheme. This in turn also implies a two approximate vertex cover maintenance 
scheme that takes 0(log n) expected amortized time per update. 

1 Introduction 



In the last decade, there has been considerable research in Dynamic Graph Algorithms where we want to 
maintain a data structure associated with some property (like connectivity, transitive closure or matching) 
under insertion and deletion of edges. Even for a simple property like connectivity, it took researchers 
considerable effort to design a polylog(n) update time algorithm |f2j [3j . In this work, we address fully 
dynamic maintenance of maximal matching in a graph. 

Let G = (V, E) be a graph on n vertices and m edges. A matching in G is a set of edges M. C E 
such that no two edges in M. share any vertex. A maximum matching is a matching that contains the 
largest possible number of edges. A matching is said to be a maximal matching if it cannot be strictly 
contained in any other matching. It is well known that a maximal matching guarantees a 2-approximation 
of the maximum matching. Ivkovic and Llyod H designed the first fully dynamic algorithm for maximal 
matching with 0((n + m) ) update time. In contrast, there exists a much larger body of work for 
maximum matching. 

Sankowski |j6l gave an algorithm for the maintaining maximum matching which processes each update 
in 0(n 1A95 ) time. Alberts and Henzinger [1J gave an expected 0(n) update time algorithm for maintaining 
maximum matching with respect to a restricted random model. Therefore the goal of a polylog(n) update 



time dynamic maximum matching algorithm appears to be too ambitious. In particular, even achieving a 
o(y / n) bound on the update time would imply an improvement of the longstanding 0{my/n) bound of the 
best static algorithm for maximum matching due to Micali and Vazirani [?]. So approximation appears to 
be inevitable if we wish to achieve really fast update time for maintaining matching. Recently, Onak and 
Rubinfeld [5 ] presented a randomized algorithm for maintaining a c-approximate (for some large constant c) 
matching in a dynamic graph that takes 0(log 2 n) amortized time for each edge update. This matching is not 
necessarily maximal, as a maximal matching would imply a factor two approximate maximum matching. In 
particular, they pose the following question - 

"Our approximation factors are large constants. How small can they be made with polylogarithmic 
update time ? Can they be made 2 ? Can the approximation constant be made smaller than two for maximum 
matching ?.." 

We resolve one of their central questions by presenting a fully dynamic algorithm for maximal matching 
which achieves O (log re) expected amortized time per edge insertion or deletion. Our bound also implies a 
similar result for maintaining a two approximate vertex cover. 

2 An overview 

Let M. denote the matching of the given graph at any moment. Every edge of M. is called a matched edge 
and an edge in E\A4 is called an unmatched edge. For an edge (tt, v) G A4, we define u to be the mate of v 
and v to be the mate of u. For a vertex x if there is an edge incident to it from M., then x a matched vertex; 
otherwise it is free or unmatched. 

In order to maintain a maximal matching, it suffices to ensure that there is no edge (u, v) in the graph 
such that both u and v are free with respect to the matching. From this observation, an obvious approach 
will be to maintain the information for each vertex whether it is matched or free at any stage. When an edge 
(u, v) is inserted, add (u, v) to the matching if u and v are free. For a case when an unmatched edge (u, v) is 
deleted, no action is required. Otherwise, for both u and v we search their neighborhood for any free vertex 
and update the matching accordingly. It follows that each update takes 0(1) computation time except when 
it involves deletion of a matched edge; in this case the computation time is of the order of the sum of the 
degrees of the two vertices. So this trivial algorithm is quite efficient for small degree vertices, but could be 
expensive for large degree vertices. An alternate approach to handling deletion of a matched edge is to use 
a simple randomized technique - a vertex u is matched with a randomly chosen neighbor v. Following the 
standard adversarial model, it can be observed that an expected deg(u) /2 edges incident to u will be deleted 
before deleting the matched edge (u, v). So the expected amortized cost per edge deletion for u is roughly 
q ^ deg(u^+deg(v) j ^ ^ deg(v) 3> deg(u), then this update time can be as bad as the one obtained by the 
trivial algorithm mentioned above; but if deg(u) is high, the update time is better. We combine the idea of 
choosing a random mate and the trivial algorithm suitably as follows. We introduce the notion of ownership 
of edges wherein we assign an edge to that endpoint which has higher degree. We maintain a partition of the 
set of vertices into two levels : and 1. Level consists of vertices which own few edges and we handle the 
updates in level using the trivial algorithm. The level 1 consists of vertices (and their mates) which own 
large number of edges and we use the idea of random mate to handle their updates. In particular, a vertex 
chooses a random mate from its set of owned edges which ensures that it selects a neighbor having a lower 
degree. This is the basis of our first fully dynamic algorithm which achieves expected amortized 0(y/n) 
time per update. 

A careful analysis of the 0(y/n) update time algorithm suggests that a. finer partition of vertices may 
help in achieving a better update time. This leads to our main algorithm which achieves expected amortized 
O (log re) time per update. More specifically, our algorithm maintains an invariant that can be informally 
summarized as follows. 



Each vertex tries to rise to a level higher than its current level if upon reaching that level, there are 
sufficiently large number of edges incident on it from lower levels. Once a vertex reaches a new level, it 
selects a random edge from this set and makes it matched. 

2.1 Related Work 

Onak and Rubinfeld also pursue an approach based on use of randomization to achieve efficient updates 
and maintain a partitioning of vertices into a hierarchy of O(logn) level that is along the lines of Parnas 
and Ron (8). The algorithm of Onak and Rubinfeld j5] takes a global approach in building level i of this 
hierarchy as follows. For level i, they consider the subgraph consisting of vertices Vi and their neighbors 
and argue that a random subset of these edges form a matching of size |Vi|/a with high probability for 
some constant a > 1. If the matching at level i falls below a predefined threshold, then a new matching is 
computed for the vertices at level i. The matching algorithm always tries to ensure that the matched edges at 
level i is always greater than | V, \ /a. As a consequence, a free vertex at level i may not get processed if the 
matching size at level i is above the threshold. This is the reason that the matching obtained is not maximal. 
The approximation factor a is an outcome of some probabilistic calculations using Chernoff bounds that is 
chosen to be a sufficiently large. Therefore, it is unlikely that any simple variation of this global approach 
can lead to a maximal matching. 

We also maintain a hierarchical partitioning of vertices but it is distinctly different from the scheme 
of Onak and Rubinfeld Q. As described earlier, the update algorithm in [5] may not process a vertex at 
level i if the matching at level i is above a certain threshold. This is a global approach of maintaining 
large matching at level i. On the other hand, we process a free vertex at level % as soon as it becomes free. 
Irrespective of the matching size at level i, we try to find a matched edge for this free vertex. This is a vertex 
centric approach for maintaining matching which ensures that the matching is maximal. Our algorithm 
achieves significantly better results than [5 ], i.e., a guaranteed factor 2 matching. The use of randomization 
is limited to choice of a random matching vertex and the O(logn) expected update time can be derived 
using 0(log n) purely random bits. 

2.2 Organization of the paper 

For a gentle exposition of the ideas and techniques, we first describe a fully dynamic algorithm for maximal 
matching that has 2 levels and achieves expected amortized 0(y/n) time per update. This is followed by 
our final fully dynamic algorithm which has log n levels and achieves expected amortized 0(log n) time per 
update (Theorem 14. lb . All logarithms in this paper are with base 2 unless mentioned otherwise. 

3 Fully dynamic algorithm with expected amortized 0{^/n) time per update 

The algorithm maintains a partition of the set of vertices into two levels. We shall use LEVEL (u) to denote 
the level of a vertex u. We define LEVEL (u, v) for an edge (u, v) as max(LEVEL(w), LEVEL (v)). 

We now introduce the concept of ownership of the edges. Each edge present in the graph will be owned 
by one or both of its end points as follows. If both the endpoints of an edge are at level 0, then it is owned 
by both of them. Otherwise it will be owned by exactly that endpoint which lies at higher level. If both 
the endpoints are at level 1, the tie will be broken suitably by the algorithm. As the algorithm proceeds, 
the vertices will make transition from one level to another and the ownership of edges will also change 
accordingly. Let O u denote the set of edges owned by u at any moment of time. Each vertex u £ V will 
keep the set O u in a dynamic hash table [7 ] so that each search or deletion on O u can be performed in worst 
case 0(1) time and each insertion operation can be performed in expected 0(1) time. This hash table is 



also suitably augmented with a linked list storing O u so that we can retrieve all edges of set O u in 0(\O u \) 
time. 

The algorithm maintains the following three invariants before the next update is processed. 

1. Every vertex at level 1 is matched. Every free vertex at level has all its neighbors matched. 

2. Every vertex at level owns less than ^fn edges at any moment of time. 

3. Both the endpoints of matched edges are at the same level. 

The first invariant implies that the matching M. maintained is maximal at each stage. A vertex u is said to 
be a dirty vertex at a moment if at least one of its invariants does not hold. In order to restore the invariants, 
each dirty vertex might make transition to some new level and do some processing. This processing in- 
volves owning or disowning some edges depending upon whether the level of the vertex has risen or fallen. 
Thereafter, the vertex will execute RANDOM-SETTLE or NAIVE-SETTLE to settle down at its new level. The 
pseudocode for insert and delete operation is given in Figure [Hand Figure |2] 

Handling insertion of an edge 

Let (u, v) be the edge being inserted. If either u or v are at level 1, there is no violation of any invariant. So 
the only processing that needs to be done is to assign (u, v) to O u if LEVEL (u) = 1, and to O v otherwise. 
This takes O(l) time. However, if both u and v are at level 0, then we execute HANDLING-INSERTION 
procedure which does the following (see Figure [T}. 

If u and v are free, then insertion of (u, v) has violated the first invariant for u as well as v. We restore it 
by adding (u, v) to A4. Note that the insertion of (u, v) also leads to increase of \O u \ and \ O v \ by one. We 
process that vertex out of u and v which owns larger number of edges; let u be that vertex. If \O u \ = y^n, 
then invariant 2 has got violated. We execute RANDOM-SETTLE(u); as a result, u moves to level 1 and gets 
matched to some vertex, say y, selected randomly uniformly from O u . Vertex y moves to level 1 to satisfy 
Invariant 3. If w and x were respectively the earlier mates of u and y at level 0, then the matching of u with y 
has rendered w and x free. So to restore invariant 1, we execute NAlVE-SETTLE(u>) and NAlVE-SETTLE(x). 
This finishes the processing of insertion of (u, v). Note that when u rises to level 1, \O v \ remains unchanged. 
Since all the invariants for v were satisfied before the current edge update, it follows that the second invariant 
for v still remains valid. 

Handling deletion of an edge 

Let (k, v) be an edge that is deleted. If (u, v) ^ M., all the invariants are still valid. So let us consider the 
nontrivial case when (u,v) G M. In this case, the deletion of (u, v) has made u and v free. Therefore, 
potentially the first invariant might have got violated for u and v, making them dirty. We do the following 
processing in this case. 

If edge (n, v) was at level 0, then following the deletion of (u, v), vertex u executes NAlVE-SETTLE(u), 
and then vertex (v) executes NAIVE-SETTLE^). This restores the first invariant and the vertices u and v are 
clean again. If edge (u, v) was at level 1, then u is processed using the procedure shown in Figure |2]which 
does the following (v is processed similarly). 

First, u disowns all its edges whose other endpoint is at level 1. If \O u \ is still greater than or equal to 
^Jn, then u stays at level 1 and executes RANDOM-SETTLE(u). If \O u \ is less than -^/n, u moves to level 
and executes NAlVE-SETTLE(u). Note that the transition of u from level 1 to leads to an increase in the 
number of edges owned by each of its neighbors at level 0. The second invariant for each such neighbor, say 
w, may get violated if | O w \ = -y/n, making w dirty. So we scan each neighbor of u sequentially and for each 



Procedure handling-insertion (u, v) 



1 if u and v are FREE then M <— M U {(u, v)}; 

2 if \O v \ > \O u \ then swap(u, v ); 

3 if |O u | = then 



foreach (u, w) G O u do 
delete (u, w) from C\,; 

X <r- RANDOM-SETTLE (u); 

if x / null then naive-settle(x); 

if w was previous mate of u then NAIVE-SETTLE(«;); 



Procedure random-settle (u) : Finds a random edge (u,v) from the owned edges of u and re- 
turns the previous mate of v 

1 Let (it, v) be a uniformly randomly selected edge from O u ; 

2 foreach (v, w) G O u do 

3 |_ delete (v , from 

4 if u is matched then 

5 X <— MATE(f); 

6 .M «- 7W\{(f,x)} 

e se 

7 |_ X NULL; 

8 „M <- M U {(u,f)}; 

9 LEVEL (n) <— 1; level(u) <- 1; 
10 return x; 

Procedure naive-settle {u) : Finds a free vertex adjacent to u deterministically 



l for each (u, x) G O u do 



if x is free then 

M <r- M U {(u,x)}; 
Break; 



Figure 1: Procedure for handling insertion of an edge (u, v) where level(u) = LEVEL (f) = 0. 



Procedure handling-deletion (u,v) 



1 foreach (u,w) G O u and level (w) = 1 do 

2 |^ move (u, w) from O u to O w ; 

3 if | O u | > y/n then 

X 4- RANDOM-SETTLE (u); 

if x ^ null then naive-settle^); 



se 



LEVEL (u) <- 0; 
NAIVE-SETTLE(li); 

foreach (u, w) G O u do 
if 1 0^ | > a/ti then 

X ^— RANDOM- SETTLE (w); 

if x / null then naive-settle(x); 



Figure 2: Procedure for processing u when (u, u) G A4 is deleted and level(u)=level(v)=1. 

dirty neighbor w (that is, \O w \ > \/n), we execute RANDOM-SETTLE^) to restore the second invariant. 
This finishes the processing of deletion of (u, v). 

It can be observed that, unlike insertion of an edge, the deletion of an edge may lead to creation of a 
large number of dirty vertices. This may happen if the deleted edge is a matched edge at level 1 and at least 
one of its endpoints move to level 0. 



3.1 Analysis of the algorithm 

We analyze the algorithm using the concept of epochs, which we explain as follows. While processing the 
sequence of insertions and deletions of edges, some matched edges become unmatched and some unmatched 
edges become matched. 

Definition 3.1 Consider any edge (u, v), and let it be a matched edge at time t. Then the epoch of(u, v) is 
the maximal continuous time period containing tfor which it remains in Ai. 

The entire life span of an edge (u, v) consists of a sequence of epochs of (u, v) separated by the contin- 
uous periods when (n, v ) is not matched. 

The third invariant implies that epochs can be partition into epochs at level and level 1. In order to 
bound the computation time per update, we calculate the computation involved in an epoch at level and 1. 
We will see that the major computations are done either at time of the start of an epoch (when the edge is 
included in the matching) or at the end of an epoch(when the edge removed from the matching). An update 
steps may involve starting and ending of many epoch. Let the total time taken by our algorithm at update 
step t is Ct- For the analysis we will redistribute this computational time among the various epochs at step 
t. Precisely, if an epoch of (u, v) starts at time t, let the computational cost associated with the start of this 
epoch be S^ u ' v \ Similarly, let the computational cost associated with the end of epoch of (w, y) be Q[ w ' y \ 
We redistribute Ct in such a way that the sum of all the St's and Q t 's is Ct at update step t. Now, we look 
at all the procedures of our algorithm and try to distribute their computation cost among various epochs. 

1. naive-settle(u) 

In the procedure NAlVE-SETTLE(u), either u finds a mate v or it becomes free. In both the cases the 



time required is 0(y/n) as u searches atmost yfn owned edges at level 0. Let us deal with these cases 
separately. 

(a) u manages to find a new mate v 

In this case, an epoch of (u, v) starts. The 0(y/n) time is associated with the start of the epoch 

of (u, v). 

(b) u becomes free 

Even in this case the time required is 0(y/n). We have to associate this cost to some epoch. 
Let us now look at the procedure naive-settle^) closely, naive-settle(m) can be called 
from two places: HANDLING-INSERTION (line 7 and line 8) and HANDLlNG-DELETlON(line 5 
and line 11). In these case NAlVE-SETTLE(u) is executed because of the following reason: the 
previous mate of u(say it?) has moved to level 1. w moves to level 1 either to execute RANDOM- 
SETTLE or due to the fact that some vertex at level 1 executed RANDOM-SETTLE and chose w 
as its random mate. So the computation cost of 0(y/n) is associated with the end of epoch of 
(u,w). 

2. random-settle(m) 

In RANDOM-SETTLe(m), u finds a random mate v from level 0. It then pulls v to level 1 to satisfy the 
third invariant. In this process, v becomes the sole owner of all the edges that have other endpoint at 
level 0(line 2,3). Since v was the owner of atmost y/n edges, the total computation time involved in 
performing this step is 0(y/n). Other steps in RANDOM-SETTLE can be executed in 0(1) time. So 
the total computation time is 0(y/n). This computation time is associated with the start of the epoch 
(u,v). 

3. handling-insertion(u, v) 

We have already taken care of line 6,7 and 8 of INSERTION-HANDLE as it calls either RANDOM- 
SETTLE or NAIVE-SETTLE on those lines. The only non-trivial step in this procedure left is line 4 and 
5. Here, vertex u is about to execute RANDOM-SETTLE(u). In a preparation for this, it becomes the 
sole owner of its edges at level 0. Since u was the owner of atmost yfn edges, the total computation 
time involved in performing this step is 0{y/n). This computation time is associated with the start of 
the epoch (u, v). 

4. HANDLlNG-DELETlON(u, v) Again note that we have already accounted for the major part of this 
procedure. Only line 1 and line 2 are the non-trivial section of this procedure that are not accounted. 
Here a matched edge of u was deleted and u starts this procedure by disowning all its owned edges 
whose other endpoint is at level 1. The number of such edges can be 0(n). Therefore, the computation 
time is 0(n) and it is associated to the end of the previous epoch of u. 

Excluding the updates that cause the start and end of an epoch of (u,v), every other edge update on 
u and v during the epoch is handled in just O(l) time. Therefore, we shall focus only on the amount of 
computation associated with the start and end of an epoch. Let us now analyze the actual computation 
associated with the epoch at level and level 1 . 

• Epoch at level 

The only non-trivial computation associated with the start of the epoch is in NAIVE-SETTLE. We have 
already seen that the computational cost of NAIVE-SETTLE is 0(y/n) and can be associated with the 
start or end of some epoch. 

• Epoch at level 1 

Consider the epoch of (it, v) at level 1. There are two ways in which epoch at level 1 starts: 
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Figure 3: Epochs at level and 1; the creation of an epoch at level 1 can destroy at most two epochs at level 
0. 

- In HANDLING-INSERTION 

The computation cost associated with the start of the epoch of (u, v) is atmost the computation 
cost required in executing the procedure HANDLING-INSERTION. We have already calculated 
that line 4,5 of HANDLING-INSERTION takes 0{y/n) time. All this computation cost is associ- 
ated to epoch of (u, v). 

- In HANDLING-DELETION 

In HANDLING-DELETION, a epoch of (u, v) may start due to calling RANDOM-SETTLE at line 
4 or at line 10. In both the cases the cost associated is equal to the cost of executing RANDOM- 
SETTLE which is 0(a/«)- 

Now let us calculate the computation cost associated with epoch of (u, v) when it ends. The only 
place where computation cost is associated with the end of the epoch is at line 1-2 of HANDLING- 
DELETION. This cost is atmost 0(n). 

From our analysis, it follows that the amount of computation associated with an epoch at level is 
0{y/n) and to level 1 is 0(n + y/n) which is 0(n). 

An epoch corresponding to some edge, say (u, v), ends because of exactly one of the following causes. 
(i) if (u, v) is deleted from the graph. 

(ii) u (or v) get matched to some other vertex leaving its current mate free. 

An epoch will be called a natural epoch if it ends due to cause (i); otherwise it will be called an induced 
epoch. Induced epoch can end prematurely since, unlike natural epoch, the matched edge is not actually 
deleted from the graph when an induced epoch ends. 

It follows from the algorithm described above that every epoch at level 1 is a natural epoch whereas 
an epoch at level can be natural or induced depending on the cause of its termination. Furthermore, 
each induced epoch at level can be associated with a natural epoch at level 1 whose creation led to the 
termination of the former. In fact, there can be at most two induced epochs at level which can be associated 
with an epoch at level 1. It can be explained as follows (see Figure [3]). 

Consider an epoch at level 1 associated with an edge, say (u, v). Suppose it was created by vertex u. If 
u was already matched at level 0, let w ^ v be its mate. Similarly, if v was also matched already, let x ^ u 
be its current mate at level 0. So matching u to v terminates the epoch of (u, w) as well as the epoch of edge 
(v, x) at level 0. We charge the overall cost of these two epochs to the epoch of (u, v). We have seen that 
the computational cost associated with an epoch at level is 0{^/n). So the overall computation charged to 
an epoch of (u, v) at level 1 is 0(n + which is 0(n). 



Lemma 3.1 The computation charged to a natural epoch at level 1 is 0{n) and the computation charged 
to a natural epoch at level is 0(y/n). 

In order to analyze our algorithm, we just need to get a bound on the computation charged to all natural 
epochs at level and level 1 during a sequence of updates. In particular, we need to bound the computation 
charged to all the natural epochs which either end during the updates or remain alive at the end of all the 
updates. By alive, we mean that matched edges corresponding to these epoch remain in the matching after 
all the update ends. 



3.1.1 Bounding the computation charged to the natural epochs destroyed 

Let t be the total number of updates. Each natural epoch at level which is destroyed can be assigned 
uniquely to the deletion of its matched edge. Hence it follows from Lemma l3TT1 that the computation charged 
to all natural epochs destroyed at level during t updates is 0(ty/n). 

Now we shall analyze the number of epochs destroyed at level 1 . Consider an epoch at level 1 initiated 
by some vertex, u by selecting a random edge (u, v) from O 1 ™ 1 . From the invariant, 0™ lt has at least y/ri 
edges. Consider the sequence of edge deletions following the creation of the epoch initiated by u. Let t u 
denote the subsequence of edges from 0" nt , and observe that the this epoch is completely determined by t u 
and the epoch ends when (u, v) is deleted. Since the vertex u selected a matched edge out of 0" nt uniformly 
at random, the matched edge of u may appear anywhere in this subsequence with equal probability 0. The 
duration of the epoch initiated by u is the number of edges in 0™ u deleted during the epoch (before (tt, v)) 
and is a random variable uniformly distributed in the range [1, See Figure H] 

Henceforth, t u will denote the truncated subsequence and we denote the duration of this epoch by \t u \. 
Since an edge e can be deleted and inserted multiple times in an update sequence, let e 1 denote the i-th 
occurrence of the edge e - e J+1 can appear only after e l is deleted. Let t u denote the labelled version of 
t u where the edges carry the additional information of their occurrence index. Since an edge can be owned 
exactly by one initiator at any time, it implies the following observation. 

Lemma 3.2 t w and t u are disjoint for w ^ u. 

In our subsequent discussion, we will not distinguish between t u and r n and assume that i u 's are disjoint. 

sequence of deletion of edges incident on u 











pi 










6 ; 







/ninit 



matched edge of u 




(i) 



sequence of deletion of edges 0™ %t 
(ii) 



Figure 4: (i) Edges incident on u and the set 0™ lt at any moment of time; (ii) The matched edge can appear 
any where in the subsequence of deletions of 0™ lt . 

For a vertex u, let t' u denote the concatenation of successive t u for all the epochs initiated by u. For a 
deletion sequence t' u , we will find the expected number of epochs which created this deletion sequence. We 



'assuming that all the edges appear and only the first deletion of any edge is relevant 



will prove that on expectation, the number of epochs for a deletion sequence t' u is 0(\t' u \/y/n). Since the t' u 
are disjoint (Lemma [3T2T >. the expected total number of epochs is 0{t/y/n), where t is the total number of 
edges deleted at level 1. From Lemma l3Tl the expected total cost charged to these epochs is 0(ty/n). 

Let X t be the random variable denoting the number of epochs at level 1 destroyed during a sequence of 
t updates. As described above X t = J2 U ^-W |> where t' u is the deletion sequence for u as described above. 
Consider the first epoch of u to be the first epoch to terminate during a sequence of t' u updates such that 
\t' u \ = k. At the moment of its creation, u owned 0* t m< edges and the duration of this epoch is a random 
variable Z distributed uniformly in the range [1...0 1 ™ 1 ]. Using the law of conditional expectation, that is, 
E[X k = E[E[X k \Z]}, and setting y k = E[X k ] and M = \0™ u \, we obtain the following recurrence. 

( £E£i(!/k-i + l) forfc>M 
Vk = { (1) 
I sEt1(!/. + l) for2<A;<M 

where y\ = Pi[Z = 1] = 1/M. For k < M, by subtracting the recurrence of y k -i from y k , we obtain 
y k — y k -i = jj (Vk-i + !)• Solving the recurrence, we obtain 




~ M 

We will prove that for for all k, y k < bek/M. We prove this by using induction on k. The base cases 
is true : VO < i < M, yi < 5ei/M . By induction hypothesis, Vi < k — 1, yi < hei/M . To show that the 
claim holds for y k , we substitute the values of y k -i, Vk-2, • • • > Dk-M i n Equation 1. We obtain 

Vk < ^(E£i^ + l) 

= l + 5efc/M-rftE£i<) 
= l + 5efc/M-(^MM±i)) 

< 1 + 5efc/M - 1 = 5ek/M 
So E[X t ] = E u E[X lt/ J < E„5e|f u |/v^ = 5ei/v^. 
Thus we conclude the following lemma: 

Lemma 3.3 For any t > 0, E[X t ] = O (t/^/n). 

Notice that this proof does not rely on the independence of the random numbers since it is obtained by 
summing up expectation of random variables representing number of epochs destroyed over subsequences 
of updates. For the base case of a sequence of length 1, we use a random number in the range [1..M], where 
M = | O" 1 ™ 1 1 < n. To choose a random number Rm in the range [1 . .M] , we first generate a random number 
R in the range [l..n] and if R < n — (n mod M) then set Rm = [R mod M) + 1 else generate a new 



R and repeat the procedure. The reason for choosing R in this way is to ensure uniform distribution in the 
range [1..M]. For n > 2M, the expected number of repetitions is less than 2. 

If we allow total independence, then we can obtain this bound with high probability as well. Since the 
matched edge of an epoch is selected independent of other epochs, it follows that duration of each natural 
epoch at level 1 is independent of other natural epochs. Hence we can state the following lemma. 



Lemma 3.4 The probability that a given epoch at level 1 has duration at most i is bounded by ^ 
Lemma 3.5 For any given q > 1, 

Pr[X t =q}< 



( Aet V /2 



\qVnJ 



Proof: If there are q epochs destroyed during t updates, at least half of them have duration < 2t/q. Hence, 
Pr [X t = q] is bounded by the probability that there are at least q/2 epochs of duration at most ^. So, using 
Lemma [3741 and exploiting the independence among the epochs, 



q ) (jl_y /2 < ( Aet y /2 

q/2) \q^J ~ [qV^J 



^^UA^J A 

For the last inequality we used (^ 2 ) < {j/j^j = (2e) 1 . □ 
We use Lemma 1331 to analyze E[X t ] and deviation of X t . 

Lemma 3.6 For any t > X t is 0{t/\/n + log n) with very high probability. 

Proof: 

We choose qo = 4 (log n + 4et/ ' yfn). It follows from Lemma [3751 that for any q > qo, Pr [X t = q] is of 
the form b q where base b < 1/2. Hence Pr[X t > qo] is bounded by a geometric series with the first term 
< 2~ qo and the common ratio less than 1/2. Furthermore go > 41ogn, hence PrpQ > qo] is bounded by 
2/re 4 . Hence X t is bounded by Oitl^fn + log n) with high probability. □ 
Notice that the proofs of I3.6l relv heavily on the total independence of random numbers used for selecting 
random mates. Using standard techniques based on generalized Chebyshev's inequality, we can obtain the 
same bound using O(logn) way independence ([?]). 

Now, recall from Lemma l3"Tl that each natural epoch destroyed at level 1 has 0(n) computation charged 
to it. So Lemmas l3.3l and l3~6l when combined together, imply the following result. 

Lemma 3.7 The computation cost charged to all the natural epochs which get destroyed during any se- 
quence oft updates is 0(ty/n) in expectation and 0(t^/n + n log n) with high probability. 

Let us now analyze the cost charged to all those epochs which are alive at the end of t updates. Consider 
an epoch (say epoch of (n, v)) such that (u, v) is still in the matching after t updates. Note that u(as well 
as v) is involved in at most one live epoch. Since epoch of (u, v) doesn't end, the cost charged to this epoch is 



1. The computation cost of starting the epoch is equal to the edges owned by u and v at the start of epoch 
(u,v). Since each vertex is involved in at most one live epoch, this cost is at most |0* m *|, 

«:epoch of v is 
live after t updates 

where 1 0* m * | is the number of edges owned by v when the live epoch of v stalled. Note that the num- 
ber of edges owned by v is less than the number of edges added to v during t updates. If t v is the the 
total number of edges added to v, then J2 V ®v < Ylv ^ n i s i m P nes that the cost associated to all 
live epochs < ^2 V t v = 2t 



2. The computation cost associated with the induced epoch which may have ended due to epoch of 
(u, v). The cost associated with these induced epoch is charged to natural epoch of (u, v). We know 
that can be atmost two such epoch( say epoch of (u,w) and (v,y)). By Lemma [3TT1 this cost is 
0(y/n). Since the total number of updates at t, the total live epoch at after t updates can atmost be t. 
So this computation charged to live epochs due to these induced epoch is at most 0(ty/n). 

It thus follows that the computation cost charged to all the live epochs at after t updates is 0(ty/n). During 
any sequence of t updates, the total number of epochs created is equal to the number of epochs destroyed 
and the number of epochs that are alive at the end of t updates. Hence using Lemma 13.71 we can state the 
following theorem. 

Theorem 3.1 Starting with a graph on re vertices and no edges, we can maintain maximal matching for any 
sequence oft updates in 0(ty/n) time in expectation and 0(ty/n + nlogn) with high probability. 

3.2 On improving the update time beyond 0(^/n) 

In order to extend our 2-LEVEL algorithm for getting a better update time, it is worth exploring the reason 
underlying 0(\fn) update time guaranteed by our 2-LEVEL algorithm. For this purpose, let us examine the 
second invariant more carefully. Let a(n) be the threshold for the maximum number of edges that a vertex 
at level can own. Consider an epoch at level 1 associated with some edge, say (u, v). The computation 
associated with this epoch is of the order of the number of edges u and v own which can be 0(n) in the 
worst case. However, the expected duration of the epoch is of the order of the minimum number of edges u 
can own at the time of its creation, i.e., 0(a(n)). Therefore, the expected amortized computation per edge 
deletion for an epoch at level 1 is 0(n/a(n)). Balancing this with the a(n) update time at level 0, yields 
a(n) = y/n. 

In order to improve the running time of our algorithm, we need to decrease the ratio between the max- 
imum and the minimum number of edges a vertex can own during an epoch at any level. It is this ratio 
that actually bounds the expected amortized time of an epoch. This insight motivates us for having a finer 
partition of vertices : the number of levels should be increased to 0(log re) instead of just 2. When a vertex 
creates an epoch at level i, it will own at least 2 l edges, and during the epoch it will be allowed to own at 
most 2 l+1 — 1 edges. As soon as it starts owning 2* +1 edges, it should migrate to higher level. Notice that 
the ratio of maximum to minimum edges owned by a vertex during an epoch gets reduced from \fri to a 
constant. 

We pursue the approach sketched above combined with some additional techniques in the following 
section. This leads to a fully dynamic algorithm for maximal matching which achieves expected amortized 
0(log re) update time per edge insertion or deletion. 

4 Fully dynamic algorithm with expected amortized O(logn) time per up- 
date 

Like the 2-LEVEL algorithm, we maintain a partition of vertices among various levels. We iterate the differ- 
ence in the partition vis-a-vis 2-LEVEL algorithm. 

1. The fully dynamic algorithm maintains a partition of vertices among [log raj + 2 levels. The levels 
are numbered from — 1 to Lo = [l°g n \ ■ We will see that the algorithm moves a vertex to level i if the 
vertex is the owner of atleast 2 l edges at that moment. So a vantage point is needed for a vertex that 
does not own any edge. As a result, we introduce a level - 1 that contains all the vertices that do not 
own any edge. 



2. We use the notion of ownership of edges which is slightly different from the one used in 2-LEVEL 
algorithm. In the 2-LEVEL algorithm, at level 0, both the endpoints of the edge are the owner of the 
edge. Here, at every level, each edge is owned by exactly one of its endpoints. In particular, the 
endpoint at the higher level always owns the edge. If the two endpoints are at the same level, then the 
tie is broken appropriately by the algorithm. 

Like the 2-LEVEL algorithm, each vertex u will maintain a dynamic hash table storing the edges O u 
owned by it. In addition, the generalized fully dynamic algorithm will maintain the following data structure 
for each vertex u. For each i > LEVEL (u), let E % u be the set of all those edges incident on u from vertices 
at level i and are not owned by u. For each vertex u and level i > LEVEL (u), the set £ % u will be maintained 
in a dynamic hash table. However the onus of maintaining 6 % u will not be on u. For any edge (tt, v) G £ % u , v 
will be responsible for the maintenance of (u, v) in £ * since (u, v) G O v . 

4.1 Invariants and a basic subroutine used by the algorithm 

As can be seen from the 2-level algorithm, it is more cost-effective for each vertex u to get settled at a higher 
level once it owns a large number of edges. Pushing this idea still further, our fully dynamic algorithm will 
allow a vertex to rise to a higher level if it can own sufficiently large number of edges after moving there. In 
order to formally define this approach, we introduce an important notation here. 
For a vertex v with LEVEL (v) = i, 



In other words, for any vertex v at level i and any j > i, 4> v {j) denote the number of edges which v can 
own if v rises to level j. Our algorithm will be based on the following strategy. If a vertex v has (j) v (j) > 2 J , 
then v would rise to the level j. In case, there are multiple levels to which v can rise, v will rise to the highest 
such level. With this key idea, we now describe the three invariants which our algorithm will maintain. 

1 . Every vertex at level > is matched and every vertex at level - 1 is free. 

2. For each vertex v and for all j > level(u), 4>v{j) < 2 J holds true. 

3. Both the endpoints of the matched edge are at the same level. 

By definition, vertices at level -1 cannot be owner of any edge. We will see that our algorithm moves a 
vertex from level -1 to a higher level as soon as it becomes the owner of some edges. In fact, the second 
invariant bounds us to act in such a manner. If a vertex v still remains at level -1 after owning some edges 
then 0„(O) > 1 > 2°, violating the second invariant. This, together with the fact that an edge is owned by 
one of its endpoints, implies that a vertex at level -1 cannot be a neighbor of another vertex at that level. 
This fact together with the first invariant imply that the matching maintained by the algorithm will indeed 
be a maximal matching. Figure |5]depicts a snapshot of the algorithm. The second invariant captures the key 
idea described above - after processing every update there is no vertex which fulfills the criteria of rising. 
An edge update may lead to violation of the invariants mentioned above and the algorithm basically restores 
these invariants. This may involve rise or fall of vertices between levels. Notice that the second invariant 
of a vertex is influenced by the rise and fall of its neighbors. We now state and prove two lemmas which 
quantify this influence more precisely. 




Lemma 4.1 Consider a vertex ufor which both invariants hold. The rise of any neighbor, say v, cannot 
violate the second invariant for u. 




-1 



matched edge 



unmatched edge 



Figure 5: A snapshot of the algorithm on Kg: all vertices are matched( thick edges) except vertex x at level 
-1. Vertex v is the owner of just the edge (v, x). 4> v (2) = 2 < 2 2 and 4> v (3) = 4 < 2 3 , so v cannot rise to a 
higher level. 

Proof: Let LEVEL (u) = k. Since the invariants hold true for u before rise of v, so 4> u (i) < 2 l for all i > k. 
It suffices if we can show that <p u (i) does not increase for any i due to the rise of v. We show this as follows. 

Let vertex v rises from level j to £. If £ < k, the edge (u, v) continues to remain in O u , and so there is 
no change in cj) u (i) for any i. Let us consider the case when £ > k. The rise of v from j to £ causes removal 
of (u, v) from O u (or Si if j > k) and insertion to £„. As a result 4> u (i) decreases by one for each i in 
[max(j, k) + 1, £], and remains unchanged for all other values of i. □ 

Lemma 4.2 Consider a vertex ufor which both invariants hold. Then any fall of one of its neighbors, say 
v, from level j to j — 1 increases 4>u(j) by at most one. 

Proof: Let level(u) = k. In case k > j, there is no change in <fi u (i) for any i due to fall of v. So let us 
consider the case j > k. In this case, the fall of v from level j to j — 1 leads to the insertion of (u, v) in 
Ei~ l and deletion from £h. Consequently, <p u (i) increases by one only for i = j and remains unchanged for 
all other values of i. □ 

In order to detect any violation of the second invariant for a vertex v due to rise or fall of its neighbors, 
we shall maintain {(j) v (i)\i < Lo} in an array 4> v \\ of size Lo + 2. The updates on this data structure during 
the algorithm will involve the following two types of operations. 

• DECREMENT-i^>(u, I): As per Lemma l4~Tl when a neighbor of v rises from level j to I, 4> v (i) decreases 
by one for all i in interval [max(j, level(u)) + 1, £]. This operation decrements (f> v (i) by one for all 
i in interval I = [max(j, level(v)) + 1, £]. 

• INCREMENT- (f>(v,i): this operation increases 4> v {i) by one. This operation will be executed when 
some neighbor of v falls from i to i — 1. 

It can be seen that a single DECREMENT-^(u, J) operation takes 0(|/|) time which is O(logn) in the worst 
case. On the other hand any single INCREMENT- ^(u , i) operation takes O(l) time. However, since <j) v (i) is 
initially and is non-negative always, we can conclude the following. 



Lemma 4.3 The computation cost of all DECREMENT-^) operations over all vertices is upper-bounded by 
the computation cost of all INCREMENT-^) operations over all vertices during the algorithm. 

Observation 4.1 It follows from Lemma \4.3\ that we just need to analyze the computation involving all 
INCREMENT-^) operations since the computation involved in DECREMENT-</>() operations is subsumed by 
the former. 



Procedure generic-random-settle (u,i) 



1 if level (it) < i then 



I lu owns edges till it reaches level i 



for each j = LEVEL (it) to i — 1 do 
for each (it, w) G £„ do 

transfer (it, w) from £i to E l w \ 
transfer (it, id) from O w to O u ; 
DECREMENT-^(u), [j + l,i]); 



7 Let (it, v) be a uniformly randomly selected edge from O u ; 

8 if v is matched then 



y 
10 



X MATE(l7); 

M<-M\{(v,x)} 

NULL 



se 



n |_ : 

12 for each j = LEVEL (l?) to i — 1 do / /v rises to level i and thus owns edges incident 
from vertices at levels level(li) to i — 1 



13 
14 
15 
16 



for each (v, w) £ £i do 
transfer (u, w) from to 
transfer (v, w) from 0„, to O v ; 
DECREMENT-(/)(l(;, [j + l,i]); 



n M ^ ML) {(u,v)}; 

18 LEVEL (it) ^— i; LEVEL(l>) 

19 return x; 



Figure 6: Procedure used by a free vertex u to at level i. 

If any invariant of a vertex, say u, gets violated, it might rise or fall, though in some cases, it may still 
remain at the same level. However, in all these cases, eventually the vertex u will execute the procedure, 
GENERIC-RANDOM-SETTLE, shown in Figure [6] This procedure is essentially a generalized version of 
random-settle(u) which we used in the 2-level algorithm. generic-random-settle(/u, i) starts with 
moving it from its current level (LEVEL (it)) to level i. If level i is higher than the previous level of it, it 
acquires the ownership of all the edges whose endpoints lie at the level G [level (it) , i — 1] . For each 
such edge (n, v) that is now owned by u, we perform DECREMENT-<^>(i>, [level (v) + 1, i]) to reflect that 
the edge is now owned by vertex u which has moved to level i. Henceforth, the procedure then resembles 
RANDOM-SETTLE. It finds a random edge (it, v) from O v and moves v to level i. The procedure returns the 
previous mate of v, if v was matched. 

Lemma 4.4 Consider a vertex it that executes GENERlC-RANDOM-SETTLE(it, i) and ends up selecting a 
mate v. Excluding the time spent in DECREMENT—^ operations, the computation time of this procedure is 



of the order of \ O u \ + \O v \ where O u and O v is the set of edges owned by u and v just at the end of the 
procedure. 

4.2 Handling edge updates by the fully dynamic algorithm 

Our fully dynamic algorithm will employ a generic procedure called PROCESS-FREE- VERTICES*}. The 
input to this procedure is a sequence S consisting of ordered pairs of the form (x, k) where x is a free 
vertex at level k > 0. Observe that the presence of free vertices at level > implies that matching M. 
is not necessarily maximal. In order to preserve maximality of matching, the procedure PROCESS-FREE- 
VERTICES restores the invariants of each such free vertex till S is emptied. We now describe our fully 
dynamic algorithm. 

Handling deletion of an edge 

Consider deletion of an edge, say (u,v). For each j > max(LEVEL(u), LEVEL (v)), we decrement 4> u {j) 
and (j) v (j) by one. If (u,v) is an unmatched edge, no invariant gets violated. So we only delete the edge 
(u,v) from the data structures of u and v. Otherwise, let k = level(-u) = level(-u). We execute the 
Procedure PROCESS-FREE-VERTICES(((li, k), (v, k))). 

Handling insertion of an edge 

Consider the insertion of an edge, say (u, v). We check if the second invariant has got violated for either of u 
or v. The invariant may get violated for u (likewise for v) if there is any integer i > max(LEVEL(w), level(v)), 
such that 4> u (i) was 2 l — 1 just before the insertion of edge (u, v). In case there are multiple such integers, 
let i max be the largest such integer. We increment (f) u (l) and 4> v (£) by one for each I > i max - To restore 
the invariant, u leaves its current mate, say w, and rises to level i max . We execute GENERIC-RANDOM- 
SETTLE(u, i max ), and let x be the vertex returned. Let j and k be respectively the levels of w and x. Note 
that x and w are two free vertices now. We execute process-free-vertices(((:e, k), (w,j))). 

Suppose that the insertion of edge (u, v ) violates the second invariant for both u and v. Wlog assume 
that u is at a higher level than v. We add (u, v) to O u and ^ EVEL ( U ) highest level to which u can 
rise after this insertion is j. This implies O u + S LEV EL(«)<fe<j \^u\ = j ust a f ter trns e dge is inserted. 
Similarly, v may rise atmost to level / and O v + Slevel(d)<£:<« \^v \ = 2 • If j > U we select u which can 
rise to the higher level and restore its second invariant. After moving u to level j, edge (u, v) becomes an 
element of Si. So S L EVEL(i>)<fc<i \^v \ decreases by 1. So, O v + Slevel^^akz \^u\ decreases by one and 
is now strictly less than 2 l , thus the second invariant for v is also restored. If i > j, then we first process v 
and move it to a higher level and the invariant of u is restored automatically. So we only process the vertex 
which can move to a higher level in this case. 

4.2.1 Description of Procedure process-free-vertices 

The procedure receives a sequence S of ordered pairs (x, i) such that x is a free vertex at level i. It processes 
the free vertices in the decreasing order of their levels starting from Lo- We give an overview of this 
processing at level i. For a free vertex at level i, if it owns sufficiently large number of edges, then it settles 
at level i and gets matched by selecting a random edge from the edges owned by it. Otherwise the vertex 
falls down by one level. Notice that the fall of a vertex from level i to i — 1 may lead to rise of some of 
its neighbors lying at level < i. However, as follows from Lemma 1431 for each such vertex v, only <p v (i) 
increases by one and <p v () value for all other level remains same. So the second invariant may get violated 
only for 4> v (i). This implies that v will rise only to level i. After these rising vertices move to level i(by 



Procedure process-free-vertices (S) 



for each (x,i) G S do ENQUEUE(Q[i], x); 
for i = Lo to do 
while Q[i] / 0do 

U «- DEQUEUE(Q[i]); 

if FALLiNGff ) then 
LEVEL (u) <r- i - 1; 
ENQUEUE(Q[i - l],f); 
for eac/j m G O s do 

transfer (it, v) from £^ to 
INCREMENT- (f)(u,i); 
INCREMENT- (/>(u, i); 
if </>«(*) > 2* then 

X «- GENERIC-RANDOM-SETTLE(u, i)\ 
if x / NULL then 

£ <r- level(x); 

ENQUEUE(Q[£],x); 



/ / v falls to i ■ 



/ / u rises to i 



e se 



/ lv settles at level % 



X «- GENERIC-RANDOM-SETTLE(f , i); 
if x / NULL then 
^ <(— level(x); 

ENQUEUE(Q[£],x); 



Function falling (v) 



1 i 4- level(w); 

2 for each (it, v) G 0„ swc/j LEVEL(u) = / do 

3 transfer (it, u) from O v to O u ; 

4 transfer (u, v) from £^ to £*; 

5 if |0J < 2* then return TRUE else return FALSE; 



1 1 v disowns all edges at level i 



Figure 7: Procedure for processing free vertices given as a sequence S of ordered pairs (x, i) where x is a 
free vertex at levelz. 



executing GENERIC-RANDOM-SETTLE), we move onto level i — 1 and proceed similarly. Overall, the entire 
process can be seen as a wave of free vertices falling level by level. Eventually this wave of free vertices 
reaches level -1 and fades away ensuring maximal matching. With this overview, we now describe the 
procedure in more details and its complete pseudocode is given in Figure [7] 

The procedure uses an array Q of size L + 2, where Q[i] is a pointer to a queue (initially empty) 
corresponding to level i. For each ordered pair (x,k) G S, it inserts x into queue Q[k]. The procedure 
executes a for loop from Lo down to where the ith iteration extracts and processes the vertices of queue 
Q[i] one by one as follows. Let v be a vertex extracted from Q[i]. First we execute the function FALLlNG(v) 
which does the following, v disowns all its edges whose other endpoint lies at level i. If v owns less than 2 l 
edges then it is decided that v has to fall, otherwise v will continue to stay at level i. 

1. v has to stay at level i 

v executes generic-random-settle and selects a random mate, say w, from level j < i (if w is 
present in Q[j] then it is removed from it and is raised to level i). If x was the previous mate of w, 
then x is a falling vertex. Vertex x gets added to Q[j]. This finishes the processing of v. 

1. v owns less than 2 l edges and has to fall 

In this case, v falls to level i — 1 and is inserted to Q[i — 1]. This fall leads to increase <p u (i) by one 
for each neighbor u of v lying at level lower than i (see Lemma |4~2~1 ). In case <j) u (i) has become 2\ u 
has to rise to level i and is processed as follows, u executes GENERIC-RANDOM-SETTLE and selects 
a random mate, say w from level j < i. If w was in Q[j] then it is removed from it. If x was the 
previous mate of w, then x is a falling vertex, and so it gets added to queue Q[j]. 

In case 1, v remains at the same level and w moves to the level of v. This renders x free and x is added 
to the Q[j]. We want to see if the invariant of any other vertex is violated in processing v and w. Since x is 
free, the first invariant of x is violated. So x is added to the queue at its level. The processing of v does not 
change <p u for any neighbor u of v. Furthermore, the rise of w does not lead to the violation of any invariant 
due to Lemma |4~T1 In case 2, v falls to level i — 1 and due to this some vertex rise to level i. All such rising 
vertex execute GENERIC-RANDOM-SETTLE. The second invariant is not violated for other vertices except 
these vertices. As in case 1, we see that processing these rising vertices may create some free vertices which 
are duly added to the queue at their level. However, their processing does not break second invariant to any 
other vertex. 

Thus we conclude the following lemma. 

Lemma 4.5 After ith iteration of the for loop of PROCESS-FREE- VERTICES, the free vertices are present 
only in the queues at level < i, and for all vertices not belonging to these queues the three invariants holds. 

Lemma 1431 establishes that after termination of procedure PROCESS-FREE- VERTICES, there are no free 
vertices at level > and all the invariants get restored globally. We want to mention this specially for the 
second invariant. The second invariant ensures that if a vertex can rise to a higher level, it should rise. 
We have seen that many vertex may rise from one level to a higher level. The algorithm processes these 
vertices in any arbitrary order. Since every vertex acts in a local way, when it is processed, it restores its 
second invariant if it is indeed violated at that moment. Lemma 14.51 ensures that there is a stable state for 
our algorithm, when all the vertices have their invariant restored. 

4.3 Analysis of the algorithm 

Processing the deletion or insertion of an edge (u, v) begins with decrementing or incrementing 4> u {i) and 
4> v (i) for all levels i > max(LEVEL(n), LEVEL (v)). The computation associated with this task over a 



sequence of t updates will be 0(t log n). This task may be followed by executing the procedure PROCESS- 
FREE- VERTICES. We would like to mention an important point here. Along with other processing, the 
execution of this procedure involves INCREMENT-^) and DECREMENT-(/>() operations. However, as im- 
plied by Observation 14.11 the computation involving DECREMENT-0Q is subsumed by INCREMENT-^) 
operations. 

Our analysis of the entire computation performed while processing a sequence of t updates is along 
similar lines to the 2-LEVEL algorithm. We visualize the entire algorithm as a sequence of creation and 
termination of various matched epochs. All we need to do is to analyze the number of epochs created and 
terminated during the algorithm and computation associated to each epoch. 

Let us analyze an epoch of a matched edge (it, v). Suppose this epoch got created by vertex v at level 
j. So v would have executed GENERIC-RANDOM-SETTLE and selected u as a random mate from level < j. 
Note that v must be owning less than 2-? +1 edges and u would be owning at most 2-? edges at that moment. 
This observation and Lemma 1441 imply that the computation involved in creation of the epoch is 0(2^). 
Once the epoch is created, any update pertaining to u or v will be performed in just 0(1) time until the 
epoch gets terminated. Let us analyze the computation performed when the epoch gets terminated. At this 
moment either one or both u and v become free vertices. If v becomes free, v executes the following task 
(see procedure PROCESS-FREE-VERTICES in Figure U} v scans all edges owned by it, which is less than 
2 J+1 , and disowns those edges incident from vertices of level j. Thereafter, if v still owns at least 2-? edges, 
it settles at level j and becomes part of a new epoch at level j. Otherwise, v keeps falling one level at a 
time. For a single fall of v from level i to i — 1, the computation performed involves the following tasks: 
scanning the edges owned by v, disowning those incident from vertices at level i, incrementing <f> w values 
for each neighbor w of v lying at level less than i, and incrementing <p v (i) by one. All this computation is 
of the order of the number of edges v owns at level i which is less than 2 l+1 . Eventually either v settles at 
some level k > and becomes part of a new epoch or it reaches level -1. The total computation performed 
by v is, therefore, of the order of J2i=k = 0(2 J ). This entire computation involving v (and u) in this 
process is associated with the the epoch of (u, v). Hence we can state the following Lemma. 

Lemma 4.6 For any % > 0, the computation associated with an epoch at level i is 0{2 l ). 

Let us now analyze the number of epochs terminated during any sequence of t updates. An epoch 
corresponding to edge (u, v) at level i could be terminated if the matched edge (u, v) gets deleted. However, 
it could be terminated by any of the following reasons also. 

• u or v get selected as a random mate by one of their neighbors present at LEVEL > i. 

• u or its mate starts owning 2' +1 or more edges. 

Each of the above factors render the epoch to be an induced epoch. We shall assign the cost of each 
induced epoch to the epoch which led to the destruction of the former. To this objective, we now introduce 
the notion of computation charged to an epoch at any level i. Note that no epoch is created at level -1 as the 
vertices at level - 1 are always free. If i = 0, the computation charged to the epoch is the actual computation 
performed during the epoch which is O(l). For any level i > 0, the creation of an epoch causes destruction 
of at most two epochs at levels < i. It can be explained as follows: Consider an epoch at level i associated 
with an edge, say (u, v). Suppose it was created by vertex u. If u was already matched at level j(j < i), 
let w ^ v be its mate. Similarly, if v was also matched already, let x ^ u be its current mate at level k. 
So matching u to v terminates the epoch of (u, w) and (v, x) at level j and k respectively. We charge the 
overall cost of these two epochs to the epoch of (u, v). 

The computation charged to an epoch at level i > is defined recursively as the actual computation cost 
of the epoch and the computation charged to at most two epochs destroyed by it at level < i. Let Q be the 
computation charged to an epoch at level i. The epoch terminated by this epoch can be at level i — 1 in the 



worst case. Also the computational cost associated with an epoch at level i is c.2 l where c is a constant. So 
we get the following recurrence: C\ < 2Cj_i + c.2\ This implies C* = 0(i2 l ). 

Lemma 4.7 The computation charged to a natural epoch at level i is 0(i2 % ). 

Henceforth we just proceed along the lines of the analysis of our 2-LEVEL algorithm analogous to the proof 
of Lemma |3.3i Let X t (i) be the random variable denoting the number of natural epochs at level i terminated 
during a sequence of ti updates. Recall that X ti = Y2 U ^|t«(i)|' wnere t u (i) is the deletion sequence for u 
By substituting M = 2*, in the the earlier analysis for two levels, using identical arguments, it follows that 

Note that Lemma 1331 shows that Pr[X t = q] < (|ffj . where M = \fn. Similarly by using 

M = 2\ we can show that if ti is the deletion sequence at level i and X tl is the random variable denoting 

the total number of epoch ending at level i, then P[X tl = q] < \^tj ■ Similar to Lemma [3T6l we 

choose qo = 4 (logn + 46^/2*). For any q > qo, Pr[X ti = q] is of the form b q where base b < 1/2. 
Hence Pr[X tl > qo] is bounded by a geometric series with the first term < 2~ qo and the common ratio less 
than 1/2. Furthermore qo > 41ogn, hence Pr[X ti > qo] is bounded by 2/n 4 . Hence X ti is bounded by 
0(t/2 l + log n) with high probability. 

Lemma 4.8 The expected number of natural epochs terminated at level i is 0(ij/2*) and 0(ti/2 l + log n) 
with high probability. 

It thus follows from Lemma 14.71 and Lemma 14.81 that the computation charged to all natural epochs ter- 
minated at level i is 0(iti) in expectation. Summing up for all levels, the expected total number of nat- 

log n 

ural epoch = it{ < logn^^"^ = O(tlogn). Similarly the total number of epochs at level i is 

i=l 

0(iti + i2*logn) with high probability. Summing up for all the levels, and using the union bound, the 

log n 

total number of epochs over all the log n levels = (it ; + i2 l log n) < t log n + log 2 n Y^i=\ ^ ^ 

i=l 

0(t log n + n log 2 n) with high probability. We can summarize as follows. 

Lemma 4.9 For any sequence of t updates, the computation charged to all the natural epochs which get 
terminated is 0(t log n) in expectation and 0(t log n + n log 2 n) with high probability. 

Let us now analyze the cost charged to all those epochs which are alive at the end of t updates. Consider 
an epoch (say epoch of (u, v) at level i) such that (u, v) is still in the matching after t updates. 

By Lemma 14771 computation charged to epoch of (n, v ) is 0(i2 l ). If 0™* was the edges owned by v 
at the start of the epoch, then since 1 0* m * | > 2* at the start of the epoch. We can say that the computation 
charged is 0(i\0™ lt \) < 0((log n)t v ) where t v is the total updates on v. Since a vertex can be part of only 
one live epoch, the total computation cost charged to all the live epochs is Y2 V 0(t v log n) = 0(t log n). 

Hence we can conclude the following result. 

Theorem 4.1 Starting from an empty graph on n vertices, a maximal matching in the graph can be main- 
tained over any sequence oft insertion and deletion of edges in 0(t log n) time in expectation and 0(t log n+ 
n log 2 n) time with high probability. 



5 Conclusion 



We presented a fully dynamic algorithm for maximal matching which achieves expected amortized 0(log n) 
time per edge insertion or deletion. Maximal matching is also 2-approximation of maximum matching. 

Our experiments show that for most of the inputs, the matching maintained is very close to the maximum 
matching. But it is not hard to come up with update sequence such that at the end of the sequence, the 
matching obtained is strictly half the size of maximum matching. We present one such example. Let 
G(V U U, E) be a graph such that V = {v\,V2, ... , v n } and U = {ui, U2, • • • , u n }. Consider the following 
update sequence. In the first phase, add edges between two points (t>j, vj) if there is no edge between them. 
Eventually this process ends when there is a complete graph on vertices of V. In the second phase, add 
edge (vi,Ui) for all i. Note that the degree of each m is one at the end of the updates. Let us now find the 
matching which our algorithm maintains. After the first phase of update ends, we have a complete graph on 
V. At that moment, we claim that the size of matching obtained by our algorithm is n. Indeed, size of any 
maximal matching on a complete graph of size n is n. Let (u, , Vj ) be an edge in the matching after phase 1 . 
Note that both these endpoints are at a level greater than -LA vertex in U is at level -1 as it does not have 
any adjacent edges after phase 1. When an edge is added, since Uj is at a higher level than ui, V{ 

becomes the owner of this edge. In the worst case, the second invariant of V{ is not violated after this edge 
deletion and nothing happens at this update step and U{ still remains at level -1. Using same reasoning, we 
can show that uj also remains at level -1 after the addition of edge (vj,Uj). Since both Uj and uj are free 
after these edge updates, there is a 3-augmenting path passing through edge (vi, Vj). In general, there will 
be a vertex disjoint 3-augmenting passing through all the matched vertices. This implies that the maximum 
matching is 2n after phase 2. But size of matching of our algorithm is n at that point. 

It would be a challenging problem to see if c-approximate maximum matching for c < 2 can also be 
maintained in 0(log n) update time. In particular, for maintaining a 3/2-approximate matching, we have to 
take care of all the 3-augmenting path. It is not clear how to extend our algorithm to handle 3-augmenting 
paths. 
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